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We consider a stationary fluid queue with fractional Brownian 
motion input. Conditional on the workload at time zero being greater 
than a large value b, we provide the limiting distribution for the 
amount of time that the workload process spends above level b over 
the busy cycle straddling the origin, as b — » 00. Our results can be 
interpreted as showing that long delays occur in large clumps of size 
of order b 2 ~ 1/IH . The conditional limit result involves a finer scaling 
of the queueing process than fluid analysis, thereby departing from 
previous related literature. 

1. Introduction. In the past ten years, there has been great interest in 
analyzing the performance of queues when the incoming traffic exhibits long- 
range dependence and self-similarity. This trend was motivated by measure- 
ments and statistical analysis of traffic in communication networks: after the 
initial findings of Leland et al. [23] in Ethernet traffic, numerous studies have 
shown that long-range dependent traffic is ubiquitous in high-speed commu- 
nication networks, and have offered partial explanations for the origin of this 
phenomenon (see e.g., Crovella and Bestavros [7]). 

One model that has received significant attention is that of a fluid queue 
that receives fractional Brownian motion (fBM) input — the so-called frac- 
tional Brownian storage (Norros [30]). As a traffic model, fBM is attractive 
because it is a stylized model (i.e., a low parameter tractable model) that 
is widely believed to be representative of long-range dependent, light-tailed 
(LRD-LT) traffic; see e.g., Norros [31] and Erramilli, Narayan and Willinger 
[11]. The tractability comes from the Gaussian self-similar characteristics of 
fBM, while its ability to approximate LRD-LT traffic is supported by vari- 
ous limit theorems, the majority of which consider superpositions of on-off 
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fluid sources with heavy-tailed connection length distributions; see Konstan- 
topoulos and Lin [19], Willinger et al. [39], Taqqu, Willinger and Sherman 
[35], Heath, Resnick and Samorodnitsky [14], Mikosch et al. [27] and Whitt 
[37, 38]. (Possible explanations for the origin of heavy-tailed distributions in 
communication networks are given by Fiorini, Sheahan and Lipsky [12] and 
Jelenkovic and Tan [17].) In this paper, we focus on an infinite-buffer fluid 
queue fed by fBM input. 

We are interested in the "local" behavior of the queue in the neighborhood 
of a "typical" time at which a customer (packet) experiences a long delay. 
We would like to answer questions like the following: What is the conditional 
behavior of the traffic at such times? How many customers (packets) will ex- 
perience large delays before congestion dissipates? This last question, about 
the size of the typical clump of long-delayed packets, provides a comple- 
ment to the commonly used performance measure given by the steady-state 
probability of a packet experiencing long delays: the latter does not dis- 
tinguish between having an occasional group of 10 packets suffering a long 
delay after an interval of 10 7 packets served, versus having a group of 10 
close-by packets with long delays, out of every 10 10 served. These situations 
can lead to drastically different "Quality of Service" (QoS) for some delay- 
sensitive applications (like audio traffic). The issue of dumpiness is also 
critical to the QoS in the finite-buffer counterpart of this model: applica- 
tions like video traffic are robust with respect to the loss of isolated packets, 
but very sensitive to the loss of a large clump of nearby packets. Note that 
under first-come-first-served discipline, a customer (packet) experiencing a 
large delay is equivalent to the workload exceeding a large threshold. Hence, 
we state our results below by conditioning on a large exceedance for the 
workload process. 

There are several results in the literature which are relevant to our work. 
Approximations for the steady-state exceedance probability for a queue fed 
by fBM were developed in a number of incremental steps: first a bound 
by Norros [30], then logarithmic-asymptotics by Duffield and O'Connell [9], 
sharper bounds by Massoulie and Simonian [26] and finally Narayan [28] 
obtained the exact asymptotics. The latter can also be obtained as a special 
case of a result by Hiisler and Piterbarg [15], who present exact asymptotics 
for the tail of the all-time maximum for a class of Gaussian processes with 
negative drift (note the steady-state workload is equal in distribution to the 
all-time maximum of negative drift fBM). The asymptotic behavior of the 
maximum workload over [0, t] as t — > oo was derived by Zeevi and Glynn [40]. 
The conditional path (in fluid scale) that leads to a long busy-period was 
studied by Norros [32] and Mandjes, Mannersalo and Norros [24]. O'Connell 
and Procissi [34] and Chang, Yao and Zajic [4, 5, 6] give the conditional 
path (in fluid scale) that leads to an exceedance of a large level b. The 
typical fluid path turns out to be nonlinear; this is in contrast with the case 
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of traffic that is short-range dependent and light-tailed (SRD-LT), where 
paths to a high exceedance are linear in great generality. More important 
for our purposes here, the typical fluid path to an exceedance of level b (with 
b large) is tangent at b, and does not spend any time above b. It follows that 
fluid analysis does not yield sufficient information about the distribution of 
the clump of packets that experience long delays; even the scaling behavior 
of the clump size is lost: the fluid limit only shows that the number of long- 
delayed packets must be of smaller order than b itself. Hence, to address 
this problem, we need to look at the workload process on a finer time and 
space scale. This is one aspect in which our work departs from most of the 
existing analyses for the fBM-driven queue. 

In the next section, we study a stationary version of the workload process, 
conditional on the workload at time being greater than b. We derive a 
conditional limit result for the deviations of the workload from level b, under 
appropriate scaling, as b — > oo. The relevant scaling consists in speeding up 
time by a factor of b 2 ~ l l H and compressing space by a factor of b 2H ~ l , 
where H is the Hurst parameter (self-similarity index) of the fBM input. The 
limiting process is a fractional Brownian motion with symmetric negative 
polynomial drift (namely — K\t\ 2H ), and started at a random level which is 
exponentially distributed; see Theorem 1. The result is proved under the 
assumption that 1/2 < H < 0.78; we believe the same result holds for all 
1/2 < H < 1, but our proof technique does not cover all this range. (Recall 
fBM has LRD increments for 1/2 < H < 1 and SRD increments for < H < 
1/2.) To the best of our knowledge, this is the first conditional limit result 
for the fBM-driven queue that uses a finer-than-fluid scale. The highly non- 
Mar kovian structure of fBM (i.e., the process becomes Markov only when 
its entire history is incorporated into the state) makes this computation 
particularly challenging. 

The conditional behavior described above is very different from the one 
observed in the case of SRD-LT traffic. Near the origin, the drift of the lim- 
iting process is close to zero; that is, conditional on a packet experiencing a 
large delay, the traffic's most natural tendency is not to restore itself immedi- 
ately to its equilibrium behavior (as is the case with SRD-LT traffic). Thus, 
many other packets are likely to also experience large delays. This suggests 
one will observe big clumps. We formalize this by deriving the conditional 
limit law for the total sojourn above level b during the busy-cycle straddling 
the origin (which is a proxy for the size of the clump of packets experiencing 
long delays). Consistent with the above scaling, the clump size scales up as 
as 5 _> see Corollary 2. Note that if one compares two fBM traffic 
sources, the one with higher Hurst parameter will exhibit asymptotically 
larger clump sizes, independent of the other parameters. Thus, in the words 
of Neidhardt and Wang [29] , for this performance measure "Hurst is always 
naughty." 
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We believe the qualitative structure of the conditional limit theory derived 
here is likely to be inherited by LRD-LT traffic in general. We also believe 
that the scaling behavior and qualitative structure is representative of queues 
fed by LRD-LT traffic. 

As mentioned earlier, in the finite-buffer counterpart of the model consid- 
ered here dumpiness (of the loss process in this case) is also an important 
issue. There is a long tradition of approximating finite-buffer loss models by 
infinite-buffer models and associated exceedance computations. In particu- 
lar, the exceedance probabilities results mentioned above are used to make 
buffer-size recommendations. In a companion paper, we consider a Brow- 
nian queue (i.e., the workload is modeled as regulated Brownian motion) 
and a GI/GI/1 with heavy-tails, and show that the qualitative structure of 
loss clumps is the same as that of the exceedance clumps. It is reasonable 
to expect a similar relation to be true here, suggesting that the loss clumps 
scale polynomially in the buffer size b in the finite-buffer model. 

From a traffic modeling standpoint, this paper contributes one building 
block towards the goal of developing a quantifiable measure of burstiness 
(which is of interest to the Internet traffic modeling community). Our view 
is that dumpiness of the loss process and burstiness of the incoming traffic 
are closely related; this paper's results therefore make rigorous the intu- 
ition that LRD-LT traffic is more bursty than SRD-LT traffic. In future 
work, we attempt related computations for stylized SRD-HT (short-range 
dependent, heavy-tailed) and LRD-HT (long-range dependent, heavy-tailed) 
traffic. These traffic types are also of interest in the network traffic model- 
ing context: there is both empirical evidence and theoretical developments 
suggesting that stable-motions or fractional stable-motions can approximate 
high-speed network traffic in situations in which the arguments in favor of 
Gaussian models like fBM fail to apply — see, for example, Konstantopou- 
los and Lin [20, 21], Mikosch et al. [27], Tsoukatos and Makowski [36], and 
Laskin et al. [22]. 

2. Model and main results. We consider a single-server fluid queue, that 
receives fractional Brownian motion as input. 

Let Z = (Z(t) : — oo < t < oo) be standard fractional Brownian motion 
(fBM) with Hurst parameter (self-similarity index) H. That is, Z is a mean- 
zero Gaussian process with stationary increments and continuous sample 
paths, started at Z(Q) = and with covariance structure given by 

EZ(t)Z(s) = (\t\ 2H + \s\ 2H -\t- s\ 2H )/2. 

The traffic process or arrival process to the queue, A = (A(t) :t G R), is 
modeled as fBM with drift, 



A(t) = Xt + 6Z(t), 
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t E R. Here A corresponds to the mean rate of the traffic, and for t > s, 
A{t) — A{s) represents the cumulative incoming traffic over the interval (s,i\. 
Let c > A denote the service rate, so that the so called netput process (or 
free process) is X = (X(t) : t £ R) given by 

X(t) = A{t) -ct = -^t + 6Z{t), 

where fi = c — A > 0. The workload (or buffer content) process W = (W(t) : 
is R) is then obtained by applying the regulator mapping to X , which in 
this case translates to 

W(t) = X(t)-MX(s), 

s<t 

t £ R. Thus constructed, W(t) represents the total amount of work present in 
the system at time t (the virtual waiting time), and corresponds to stationary 
regulated fBM. Note that W satisfies the relation 

W{t) = X(t) - X(s) + W{s)V (x(s) - inf X(u)), t>s. 

V s<u<t / 

We are interested in studying the "local" behavior of the workload process 
around a "typical" exceedance of level b. For this purpose, we define the 
process Y b given by 

Y\t)^b- 2H+1 (W(tb 2 - 1 / H )-b), 

i£R. Note that this is indeed a description of the local behavior near the 
origin: conditional on {W(0) > b}, the length of the busy-cycle straddling 
the origin scales up linearly in b, while the speed-up factor in Y b is b 2 ~ l l H . 
Hence, for any fixed t > 0, as b increases the path (Y b (s) : — t < s <t) rep- 
resents an ever smaller fraction of the busy-cycle straddling the origin. 

Our main result gives the limiting law of Y b conditional on Vv"(0) > 6, as 
b — ► co. 

Theorem 1. If 1/2 < H < (y/Tf - l)/4 then 
(1) (Y b \W(0)>b)^Y asb^oo, 

in C(— oo,oo) with the topology of uniform convergence on compact sets, 



|2H _|_ Q7(+\ + r- TO v(n\ II 7 V ( n\ ~R 1 
Kl ~ 2(1-H) v H > ' ana K2 — 6 i H iH 



where Y{t) = Y(0) - Kl \t\ 2H + 9Z(t), t G R, Y(0) JL Z, Y(0) = ± Exp(l), 



As mentioned above, we believe the result is likely to hold for all 1/2 < 
H < 1, but a portion of the argument in the proof only works for 1/2 < H < 
(VT7 — l)/4; see Step 12 in the proof of Theorem 1. 

It is interesting to compare the case H > 1/2 with the case H = 1/2, i.e., 
with the result for the Brownian queue. There the limiting law is that of a 
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Brownian motion with symmetric negative linear drift. Here, it is fractional 
Brownian motion with symmetric negative polynomial drift. The results 
are, to some degree, in agreement, since the conditional limiting law for 
the Brownian queue corresponds exactly to replacing formally H = 1/2 in 
Theorem 1. Qualitatively, however, there are significant differences. First, 
in the Brownian case the limiting law is obtained without the need to scale 
time and space. Also, in the Brownian case the drift is equal to — fi\t\. In 
particular, the traffic process on [0, oo) obeys its usual (unconditional) law. 
In contrast, the polynomial drift —Ki\t\ 2H obtained in the LRD case has a 
derivative equal to zero at the origin. This, together with the scaling in place, 
indicates that from time onwards, and for a period of order 6 2-1 /^, the 
traffic intensity will be close to 1. Thus, the traffic process does not revert 
immediately to its unconditioned dynamics (as in the Brownian case), but 
rather the traffic load remains higher than average, putting the queue in 
heavy-traffic, for a significant period. 

The fact that the limiting process Y has symmetric drift is also interesting, 
and perhaps somewhat surprising, since the typical path to level b in fluid 
scale is not symmetric: the buildup from to level b is slower than the 
"draining" back to 0; see Chang, Yao and Zajic [5]. 

Next, we examine the effect on the size of the clump of long-delayed 
packets. The busy-cycle straddling the origin is (tl,tr), where 

T R = mf{t>0:W(t) = 0}, 

t l = sup{i < : W{t) = 0}. 

Let V b denote the total sojourn above level b over (tl,tr), i.e., 

v b ^ ( TR i {Koo) {w{t))dt. 

Its conditional limit distribution is given in the next result. 
Corollary 2. If 1/2 < H < (y/Tf - l)/4, then 

POO 

(b- 2+1 l H V b \W{Q)>b)^V± l [0>oo) (Y(t))dt asb^oo, 

J —oo 

where Y is as in (1). 

Proof. Fix a > and let (p a : C{— oo, oo) — > R be given by ip a (y) = 
S-a 1 (o,oo) (y(*)) dt. Put G a = {y G C(-oo, oo) : lim^^ y(t) = lim t __ 00 y(t) = 
— oo, Leb({t G (—a,a) :y(t) = 0}) = 0}, where Leb denotes Lebesgue mea- 
sure. Note ipaiVb) Waiv) whenever y^ converges to y G G a uniformly on 
compact sets and P(Y G G a ) = 1; it then follows from Theorem 1 and the 
continuous mapping theorem that ((p a (Y b )\W(0) > b) ip a (Y) as b — > oo. 
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Since <fi a {Y) / V a.s. as a f oo, it follows that lim^oo lim^oo P(ip a (Y b ) > 
x\W(0) >b) = P(V > x), x E R. On the other hand, \P{cp a (Y b ) > x\W{0) > 
b) - P(b~ 2+1 / H V b > x\W(0) > b)\ < P(F(a,b)\W{0) > b), where F(a,b) = 
{3a<t< T R :Y b (t) > 0} U {3t l < t < -a:Y b (t) > 0} U {\t l \ < ab 2 ^ 1 / 11 } U 
{t r < ab 2 - 1 ' 11 }. But lim o ^ oo lim 6 _ +oo P( J F(a ) 6)|W r (0) > b) = 0, whence the 
result follows. □ 

We see that the total sojourn above level b scales as b 2 ~ l l H . In the Brow- 
nian case, H = 1/2, the unsealed total sojourn above b has a conditional 
limiting distribution. Hence, in a queue fed by fBM traffic with H > 1/2, 
larger delays tend to occur in larger clumps, whereas in the Brownian case 
the distribution of the typical exceedance clump of a given level is roughly 
independent of the magnitude of the said level. 

It is interesting to note that the random variable V plays a role in the 
prefactor of some asymptotic tail probabilities (EF _1 being a variant of 
Pickands' constant), and also in other conditional limit results for sojourns 
above a high level of stochastic processes: For example, it appears in The- 
orems 3.3.1 and 5.5.1 of Berman [3], which deal, respectively, with high 
sojourns of stationary Gaussian processes and stationary diffusions over a 
finite interval; in both results, the probability that the (scaled) sojourn above 
b is larger than x converges to ~E(1/V; V > x), where V represents, like V 
above, the total sojourn above zero of a Gaussian process with symmetric 
negative polynomial drift, started off at an exponentially distributed offset. 
(Those results are, however, quite different from that in Corollary 2: the 
workload process W is neither Gaussian nor a diffusion, we cannot restrict 
it to a finite time interval, the limiting distribution of the scaled sojourn is 
directly the distribution of V rather than a transformation of it, and in the 
case of W the sojourn above b scales up with b, whereas in Berman's results 
it scales down with b.) 

An important tool in the proof of Theorem 1 is the following result, which 
states that if a negative-drift fBM is conditioned on hitting level b, the hitting 
time of b will be "very near" its most likely value, t* = j^r^rm ■ 

Theorem 3. Let T b = M{t > : X(t) > b}. For all r > and e>0, 
P(\T b - t*\/b H+r > e\T b < oo) — 

as b — > oo . 

For the proof, we refer the reader to Lemma 7 in Dieker [8], who proves 
a more general result; Theorem 3 corresponds to the particular case of fBM 
and (in his notation) 5(u) = 0{u H+r ). 

Note that Theorem 3 is a refinement on the fluid behavior of T b obtained 
from large-deviations analysis, namely that (\T b — t*\/b\T b < oo) — > in 
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probability as b — > 00. It turns out that the fluid scale result is too coarse 
for our purposes. A further refinement can be found in Hiisler and Piterbarg 
[16], Theorem 1, who show that ((T& — t*)/b H \Tb < 00) converges weakly to 
a Gaussian random variable. 

3. Proofs. Before proceeding with the proof of Theorem 1, we recall the 
notion of asymptotically equivalent events. We say that two collections of 
events (E b : b > 0) and (E\ : & > 0) are asymptotically equivalent if 

P(E b AE b ) = o(P(E b )) 

as b — ► 00. In the proof, we will use the following well-known fact. 

Lemma 1 (Asmussen and Kliippelberg [2], Lemma 2.4(a)). If {E b : b > 0) 
and (E\ : b > 0) are asymptotically equivalent, then 

||p(-|^)-p(-|^)IHo 

as b^ 00, where \\ ■ \\ denotes total variation distance. 

Throughout this section, we use the following notation: for any process Z 
and cr-fields T and Q, we denote Z F {t) = E(Z(t)\T), Z T {t) = Z(t) - Z T (t), 
and Zg(t) = Z T {t)-E{Z :F {t)\g). The complexity of the argument requires a 
significant amount of additional notation, which is introduced at the points 
when first needed; a notation summary table is included at the end for ease 
of reference. 

Proof of Theorem 1. To find the conditional limit law of Y b given 
{W(0) > b}, we study separately its intercept with the ordinate axis and 
its deviations from its value at the origin; that is, we study the random 
variable Y b (0) = b- 2H+1 (W(0) - &)— the scaled "overshoot" of W over level 
b at time — and the process Y b = (Y b (t) :iel) given by 

yb^ A b -2H+l X ( b 2-l/H t) 

(Note Y b (t) = Y b (t) - Y b (0) for t £ (t l ,t r ).) We find the conditional limit 
laws of Y b and Y b (0) separately, and then show that they are asymptotically 
independent as b — > 00. 

Observe that {W(0) > b} = {T b < 00}, where T b = inf{t > 0:X(-t) < 
—b}. The event {T b < 00} depends on the whole history of X over (— 00, 0], 
and the law of Y b conditioned on such an event is hard to analyze. However, 
by Lemma 1, to prove that (Y b \T b < 00) ==>• Y it is enough to show that 
(Y b \A(b)) =^ Y, where A(b) is some appropriately chosen asymptotically 
equivalent event for which the conditional law is more tractable. 
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To begin with, one can restrict attention to the history of X over a fi- 
nite interval. Put t L = t* - b H+r and t R ■ = t* + b H+r ', for some small r > 
(specified later). By Theorem 3, it follows that the event 

Mb) = {3s G l-t R , -t L ] : X(s) = -b} 

is asymptotically equivalent to {T b < oo}. The asymptotically equivalent 
event A(b) that we construct is a subset of A)(fr)> an d encompasses the 
intuition that, given {T b < oo}, X will attain values close to —b in a neigh- 
borhood of —t*. 

If one formally considers the law of Y b conditioned on X(—t*) = —b, then 
it is an easy exercise to verify that the conditional mean and covariance 
functions of this Gaussian process converge to those of the desired limit 
process Y as b — > oo, where 

Y(t) = -Ki\t\ 2H + 0Z(t). 

The first several steps in the proof show that this is still true when condition- 
ing on TL — a((X(s) : — tjt < s < —tL)), as long as X remains close enough 
to — b and satisfies some regularity conditions; making precise how close and 
the requisite regularity conditions will lead to the appropriate definition of 
A(b). 

The remainder of the argument is broken into "steps." A high-level overview 
can be obtained by reading the first portion of each, with the statement of 
the respective intermediate results, while skipping the details that follow. 
It is because of conditions imposed in Step 11 that the proof needs the as- 
sumption H < (\/T7 — l)/4; the intermediate results presented in the other 
steps are valid for 1/2 < H < 1; see discussion in Step 12. 

Step 1. The conditional mean given the "endpoint" X{—ti). Define the 
event 

Ai(6)4{x(-i i )€[-6,-6 + 6'*]} J 
Ai(6) = Mb) n Ai{b), and K = <r{X{-t L )). In this step we show that if 

(2) 0<7<# + r and < r < (1 - H) 2 /H, 

then, given M > 0, there exists e\ > such that, for all large enough b, 

i Xiib) \E(Y b (s)\>c) +K1 \ t n <b-^i Ii(b) , 

a€(-M,M). 

To prove this, put II = b~ 2+1 / H ti and note that 

L var Z(-tL) J 
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b -2H+l 
b -2H+l 



{\{X(-t L ) - ^t L )(\s/t L \ 2H + 1 - (1 + s/t L f H ) ~ vb 2 - l ' H s} 



(-b-^t*)\s/t L \ 2H 



+ b 



-2H+1 



(3) 



{ - b -^\ l-{l + s/t L ? H )-^t* 

f fi^-(t* -t L ) + - [X(-t L ) + b + nif - t L )] 

x(\ s /t L \ 2H + l-(l + s/t L ) 2H ) 



„ I „\2H /+* /. \2H 
-Kl\S\ (t jt L ) 



(l-(l + s/t L ) 2H + 2Hs/t L 



2(1 -H) 



A-2H 



+ M f + L— [X(-t L ) + b + »b H+r ] 

x(\s/t L \ 2H + l-(l + s/t L ) 2H ). 
Since \X(-t L ) + b\ < W < b H+r on A^b), it follows that, on 

|E(y 6 ( s )|/c) + Kl | s | 2 ^| 

< Ki|s| 2H [(l + b H+r /t L ) 2H - 1] + ^ - (1 + s/i L ) 2H + 2Hs/t L \ 



+ 



K\\S 



b l-H+r 



tL 



2/i + (l + /x) 



2(1 -H) 
\s/t L \ 2H + !-(! + s/t L ) 2H 



2H 



2Hb H+r b 2 ^^ 



tL 



+ 



2(1 -H) 



s/t L 

H(2H-l)(s/t L ) 



hl-H+r I I 

+ __.U[2 j u + 2if(l + / u)] 

for any < e\ < (1 — H) 2 /H - r. 

Step 2. Conditional mean given both endpoints. Define 

g±a(X(-t L ),X(-t R )), 

A 2 (b) = {\Z K (-t R )\<V*>}, 

where 7/3 > is a constant (specified later), and A 2 (b) = A$(b) n j4i(6) H 
^42(6). In this step, we show that if 

(4) m <H 2 + (l-H) 2 (l + H)/H + r(2H - 1) 
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then, given M > 0, there exists £2 > such that, for all large enough b, 
I Mb) \E(Y\s)\g)-E(Y\s)\lC)\ <I Jn(h ,b-*\ 



A 2 (6) 



se(-M,M). 

To verify this, put D = 2b H+r and observe that, on A^{b), 



\E(Y b {s)\Q)-V(Y b {s)\1C)\ 



6b -2H + l 



< 



var Z(-t L ) var Z{-t R ) - cov(Z(-t L ), Z(-t R )) 2 
x | var Z(-t L ) cov(Z(sb 2 - 1/H ), Z{-t R )) 

- C ov(Z(-t L ),Z{sb 2 ~ l l H )) cov(Z(-t L ), Z(-t R ))\ 
9b- 2H+1 b r » 

4(1 + D/t L y» - (1 + (1 + D/t L fH - {D/t L y H y 



■\z K {-t R )\ 







s 


2H 




X 


'( 




+ 


Kf"( 






T L 





s 

Q bV3 -2H+l 



211 



+ 1 



2H 



D\ 2H 

1+ d 



27i 



■AH{H-1) 



±{D/t L y H - t L t L 

= 0(ft^3-H 2 -(l-//) 2 (l+H)//f-r(2//-l) j 

= o(6" e2 ) 

as 6^00, for any < e 2 < H 2 + (1 - H) 2 (l + H)/H + r(2H - 1) 

Step 3. Conditional mean given TL. Here, we show that the conditional 

mean of given 7i is close to its conditional mean given Q, as long as 

the values of Z over the interval [—t R , —t£[ remain close enough to their 

conditional expectation given Q. 
To be specific, define 

A 3 (b) = {\Z g {-t L - s)\ < m a x(s H b 5 ,s H ^),0< s < b H+r }, 

A 4 (b) = {\Z g (-t R + s)\< m^{s H b 5 , /-"'j.O < s < b H+r }, 

where rji > and 5 > will be specified later, and put 

A(b) = A {b) n Ai(b) n A 2 (b) n A3 (6) n A* (6). 

In this step, we show that if 

(5) 5 < (1 - #) 2 (1 + - (1 - H)(r + m) + rrj u 
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then, given M > 0, there exists £3 > such that, for all large enough b, 



I 



A{b) \E(Y b (s)\H)-E(Y b ( S )\g)\<I A{b) b- 



£3 



(6) 

s G (— M, M). 

(Note the term s H ~ m within the maximum controls the deviations of Z 
from its conditional expectation given Q in the immediate neighborhood of 
the endpoints — tR and — t^; the term s H b 5 controls the deviations of Z from 
its conditional expectation given Q away from the endpoints.) 

To prove (6), we find a representation of the LHS in which the differences 
Z(— ti, — s) — E{Z{— f £ — s)\Q) appear explicitly, inordertousetheregularity 
conditions contained in ^3(6) and ^4(6). For this purpose, it is convenient 
to define the process Z by Z{u) = Z(—ti + u) — Z{—tL), m6K, Note Z is 
standard fBM, and that in terms of Z, H = o(Z(t L ), (Z(u) : - D<u<0)) 
and Q = a(Z(t L ), Z(-D)), where D = 2b H+r as before. Also, 



Y b (u) 



6b~ 2H+l [Z(t L + b 2 ^ H u) - Z(t L ) - (v/0)b 2 - l ' H u} 



Put J 4 a((Z(u) : -D<u<0)). Then, for fixed s,t G (— M, M), and 6 
large enough that ^ > b 2_1 ^ H M, we have the representation 



(7) 



ZjtL + tf-Ws) = Z/itL+b 2 - 1 /* 

\z(t L + b 2 - l / H t)) \zJ(t L + b' 



+ 




where (Wi, W 2 , W 3 ) is i.i.d. standard Gaussian, W\ JL J, (W 2 , W 3 ) JL W and 
the (deterministic) constants (o«) are as needed to match the covariances 
on both sides. 
Hence, 

Y\s) = 6b~ 2H+1 [Z J (t L + b 2 ^/ H s ) - Z J (t L ) 
(8) + (a 21 /a u - l)(Z(t L ) - Z J {t L )) 

+ a 22 W 2 -(fi/9)b 2 ^ H s], 



(9) 



Zg(t L )-(a 21 /a n -l)Zg(t L )}. 



and it follows that 

E(Y b ( S )\H)-E(Y b (s)\g) 

= 9b~ 2H+l [z3{t L + b 2 - l / H i 

Gripenberg and Norros [13] show that, for u > 0, the conditional expec- 
tation Z^(u) can be expressed as 

r o 

g(D,u,v)dZ(v), 

D 



Z J (u) 
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-H+l/2 -H+l/2 fu w H - 1 / 2 {w+D) H - 1 / 2 



dw for 



where g(D,u,-v) 4 C H (D - v )-»+V* v -»+V* /« 
< v < -D and C# = sin(7r(i? — 1/2))/tt. Here, the integral against fBM can 
be defined both as an I? limit or as an almost sure limit by approximating 
g(D,u,-) by a sequence of simple functions; for a discussion of integration 
against fBM see for example Duncan, Hu and Pasik-Duncan [10]. Moreover, 
a similar argument to that in Norros, Valkeila and Virtamo [33], Lemma 2.2, 
shows that an "integration by parts" formula holds, namely 



u 



o 



g'(D,u,v) 



D 

Z{-D) 
D 



Z{v) + -Z(-D) 



dv 



(10) 



= lim 

e\0 



Z(v) + -Z(-D) 



D+e 

Z{-D) 

D J-D+e 

g{D,u,-D + e)[Z(-D + e 
J(-D) 



dv 



g(D, u, v) dv 

D 

g'(D,u,v 

g(D, u, v) dv + g(D, u, —e)Z(—e) 
Z(-D)] 
[g{D,u,-e)+g(D,u,-D + e)} 



D 

where g'(D,u,—v) = dg(D,u, — v)/dv and the limit is a.s. well defined by 
the Holder continuity of fBM paths (cf. Norros, Valkeila and Virtamo [33], 
Lemma 2.2). Similarly, using a sequence of simple functions to approximate 
g(D,u,-) one can show that ~E(Z^(u)\Q) has a representation given by the 
RHS of (10) with Z(-) replaced by Z G (-). It then follows that 



D 



g'(D, u, v)Zg(v) dv = lim 



e\0 J-D+e 



g'(D, u, v)Zg(v) dv. 



We use this representation of Zg (u) to show that, if (5) holds, then on 
A(b) the RHS in (9) is bounded by 0(b~ £s ), for some £3 > 0. The details are 
given in Lemma 2. 

Step 4. Convergence of the covariance function. If r < (1 — H) 2 /H, then 
given i,sGR, there exists Eq > such that, for all large enough b, 



cov (Y b (t),Y b (s)\H) - (e 2 /2)(\t\ 2H + \s\ 2H 



t 



s\ 2H )\<b~ £a . 



This is proved as Lemma 3 below, using the representation (7). 

Step 5. Convergence of finite- dimensional distributions of (Y b \A(b)) . Note 
that A(b) G 7i, and that (Y b \7i) is a Gaussian process. The convergence of 
the conditional mean and covariance functions from previous steps, together 
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with an application of the Cramer- Wold device allow us to conclude that 

(y 6 |i(j)) f ^f, 

where 

Y(t) = - Kl \t\ 2H + 6Z{t), 

t G K, and denotes convergence of finite-dimensional distributions. 

Step 6. Tightness of (Y b \A(b)). We now show that {(Y b \ A(b)) : b > 0} is 
tight. We do so by showing that, given M > 0, there exists Cm > such 
that, for all large enough b, 

(11) I A{b) E((Y b (v) - Y b {u)f\H) < I A{b) C M \u - v\ 2H 

for all u,v S [—M,M], which implies 

E«f > w - *. W )vw) . EJ ^' E(( ; t ;; ) ( -f V))2|w) < cmi, - »r , 

which in turn is a sufficient condition for tightness with respect to the topol- 
ogy of uniform convergence on compact sets; see e.g., Karatzas and Shreve 
[18], page 64. 

To verify (11), note that 

E((Y b (v) - Y b (u)) 2 \H) = [E{Y b (v)\H) - E{Y b {u)\H)} 2 

+ vm(Y b {v)-Y b {u)\H). 

But, conditional on 7i, (Y b (v) — Y b (u)) is Gaussian, and its conditional 
variance is bounded above by its (unconditional) variance, so that 

var(y 6 (^) - Y b {u)\H) < d 2 \u - v\ 2H . 

Also, in Lemma 4, we use the representations (8) and (10) to show that the 
function u \— > E(Y b {u)\TL) is differentiable and that there exists a constant 
Cm such that, on A(b), 

dE(Y b {u)\H) 



< Cm, 



du 

u G [—M,M]. In particular, on A{b) 

\E(Y b (v)\H) - E(Y b (u)\H)\ < C M \u - v\ 

for u,v £ [-M,M]. Putting C M = C 2 M + 4M 2 + 6 2 , (11), and hence, the 
tightness of {(Y b \A(b)) : b > 0}, follows. 

Step 7. Limit law of Y b . Given the tightness of {(Y b \A(b)) : b > 0}, and 
the convergence of its finite-dimensional distributions, it follows that 

(Y b \A(b))^Y, 
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where ==> denotes weak convergence in C{— 00,00) (with the topology of 
uniform convergence on compact sets). 

We postpone to Step 11 the proof that A(b) is asymptotically equivalent 
to {VK(0) > b}. Assuming this for the moment, it follows by Lemma 1 that 

(12) (Y b \W(0)>b)^Y. 
Step 8. Limit law ofY b (0). We now show that 

(y 6 (0)|^(0) >b)^Y(0) = — Exp(l). 

^2 

The proof uses the exact asymptotic for the tail of W(0) as given by 
Hiisler and Piterbarg [15], Corollary 2: 

(13) P(T b < 00) ~ Kb (1 ~ H)2/H {l - Qi^b 1 " 11 )), 

where K is a constant (which they provide explicitly, in terms of Pickands' 
constant), K3 = y/ ^2/(1 — H) and <£(•) is the c.d.f. of a standard Gaus- 
sian random variable. Since (1 — &(x)) ~ (p(x)/x as x — > 00, where <j)(x) = 
exp(— x 2 /2)/\/2tt is the standard Gaussian density, we obtain 



Hence, 



P(T b < 00) ~ (A7k 3 )6- (1 - h)(2 " 1/h) exp(-(«|/2)6 



lim P(Y b {0) >x\W(0) >b) 

= hm — — 7 : 

6^00 P(T b < 00) 

b k(1-H)(2-1/H) 



| 2(1-J^)^ 



: lim , 

Kl^b + Xb 211 - 1 )^- 11 ^ -b^ 1 -*)} 



x exp 



2 



= exp(-(«i/2) hm [(b + xb 211 - 1 )^ 1 -^ -b 2 ^}) 
= exp(-(K§/2).2(l-fl-)x) 

= exp( — K2X). 

Step 9. Asymptotic independence of Y b (0) and Y b . We now show that 
Y b (0) and Y b are asymptotically independent, in the sense that 

P(Y b G ., Y b {0) > x\W{0) >b) — ► P(Y G •) • P((1/k 2 ) Exp(l) > x) 

as b — > 00. 
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To verify this, note 

P(Y b e-,Y b (0)>x\W(0)>b) 

= P(Y b e -,W(0) >b + xb 2H ~ 1 )/P(W(0) > b) 
= P(Y b e -\W(0) >b + xb 2H - l )V(W{V) > b + xb 2H - 1 )/P{W{<d) > b) 
= P{Y b E -\W(0) > b)P(Y b (0) > x\W(0) > b), 
where 6 = 6 + xb 1 . But we have shown above that 

p(y fe (o) > z|w(o) > 6) -> p(y(o) > x). 

Hence, it is enough to show that P(Y b £ -\W(0) > 6) -> P(Y G •). Observe 
that Y b = U b where 

2H-1 



/ 6 \m-L . 



P(t,b) = t • (i](b)/b) 2 l l H and 77 is the inverse function of 6^6 + x6 2 ^ 1 . 
Note that (i) /?(-, 6) converges to 1 1 — »• t uniformly on compacts sets as 6 — > 00, 
M. UJa-i ) 2 *" 1 / 1 as 6 ^ 00, (iii) 6 - 00 <=> 6 - 00, and (iv) {(Y S | W(0) > 
6) : 6 > 60} is tight (because of (12)). From (i)-(iv), it follows that for arbi- 
trary e > 0, 

P(p(fA Y b ) > e\W(0) > 6) — > 

as 6 — > 00, where p is the usual metric on C/^^y Since (y b |VF(0) > 6) =>• 
Y" as 6 -» 00, it then follows that (t/ 6 |W(0) > 6) =>> Y as 6 — ► 00, i.e., 

(y fo |vr(o)>6)=^y, 

as desired. 

Step 10. Limit law ofY b . It follows from the previous three steps that 
((Y b (0),Y b )\W(0)>b)^(Y(0),Y) 

in R x C(-oo,oo), where y(0) = (1/k 2 ) Exp(l), y(0) X Z. Since y b (t) = 
y b (0) + y 6 (t), t € (tl,tr), and since — T£,tr=^oo as 6— > 00, Theorem 1 
then follows by the continuous mapping principle. 

It only remains to be verified that A(b) and {W(0) > 6} are asymptotically 
equivalent, as we claimed in Step 7 above. We verify this in the next step. 

Step 11. Asymptotic equivalence of A(b) and {T b < 00}. We need to show 
that P{A(b)A{T b < 00}} = o(P{T b < 00}). By Theorem 3, 

P({T fc <oo}\ylo) = o(P(T 6 <oc)), 



CONDITIONAL LIMIT THEOREMS FOR FBM 17 

so it is enough to show that 

P(A (b)\A(b))=o(P(T b <oo)). 

In what follows (and for the rest of the paper) we frequently write Aq, . . . , A4, 
omitting the explicit dependence on b. Note that 

A \ A(b) c (A n Af ) u (A 1 n A c 2 ) U (At n A$) U (Ai n A C A ) 

c A 5 u (A n A 6 ) u (Ai n a£) u (A 1 n A3 7 ) u (Ai n Af ), 

where A 5 = {X(-i L ) < -6} and A 6 = {X(-t L ) > -6 + 6^}. Thus, 
P(A \ A(6)) < P(A 5 ) + P(A n A 6 ) + P(Ai n A%) 

(14) 

+ P(A 1 nA 3 7 ) + P(A 1 nAJr). 

We show that all the terms on the RHS in (14) are o(P(T b < 00)) as 
b — > 00, as long as the parameters 7, r, 771 , 7/3 and 8 are chosen appropriately. 
We list the required conditions on the parameters below; the details of the 
arguments are relegated to Lemma 5. 

That P(As) = o(P(T b < 00)) follows from Theorem 3. 

The proof that P(AoH Aq) = o(P(T fe < 00)) makes repeated use of BorelPs 
inequality (Lemma 7) and requires that 

(15) -/ > H + r- (H -1/2)(1- H -r). 

The proof that P(Ai n A$) = o(P(T b < 00)) uses the additional condition 

(16) 7] 3 >H 2 + rH + (l-H + r)/2. 

The proof that P(Ai Pi A§) and P(A X n A^) are o(P(T b < 00)) relies on 
BorelPs inequality and the Holder continuity of fractional Brownian paths, 
and uses the condition 

(17) 8>(l-H + r)/2. 

Step 12. Final remarks. If 1/2 < H < (y/Vf — l)/4, then the parameters 
r, 7,773,5 and 771 in the definition of A(6) can be chosen so as to satisfy (2), 
(4), (5), (15)-(17). 

It is because of conditions (16) and (17) that the proof does not work for all 
1/2 < H < 1: these conditions, when combined with (4) and (5), require that 
(1 -H)/2<(1- H) 2 (l + H)/H, which in turn requires H < (\/Vf - l)/4. 
The need for (16) and (17) may be a consequence of the style of proof 
rather than necessary conditions. If one is able to prove that P(Ai n A§ ) 
and P(A 2 nA^) are o(P(T b < 00)) without imposing (16) and (17), then the 
proof of Theorem 1 would work for all H G (1/2, 1). □ 
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Lemma 2. Fix s > 0. If 5 satisfies (5), then there exists £3 > such that 
(18) I A{b) \E(Y\s)\H)-E(Y b ( S )\g)\<I A{b) b-^ 
for all large enough b. 



Proof. We analyze separately each of the terms on the RHS of equa- 
tion (9). Using that g'(D,u,v) = g'(D/u, l,v/u)/u (which is straightforward 
to verify), we note that on A(b), 



zi(t L )\ 









Jo 




fD/t L 







< bh L » 



\g'(D/tL, 1, -x)Zg(-xtL)\ dx 
rD/2t L 

/ \g'(D t L ,l,-x)x H -^dx 
Jo 



+ 



D/t L 
D/2t L 



f(D/t L ,l,-x)\(D/t L - x) H ~^ dx 



where the last step follows since, on A(b), \Zq(-v)\ <ma^{v H b 5 ,v H ~^) for 
< v < D/2 and \Zg(-v)\ < max((£> - v) H b s , (D - v) H ~^) for D/2 <v< 
D. Since 



g'(D,u,-v) = -(H - l/2)g(D,u, -v)[y- 1 - (D - v)' 1 } 



-C H v 

it follows that 

\zi(t L )\ 

< b 5 t L H C H 

"D/2t L 



B + i/* {D _ v) -H+i/2 f u y H ~ 1/2 (y + D) H ^ 



1/2 



(y + v) 



dy, 







.1/2-m ( £. 



x 



-H+1/2 



l y H-l/2 {y + D/tL) H-l/2 



{y + xf 



dydx 



[ D/tl 

JD/2t L 



-H+1/2 



x 



l yH-l/2 {y + D/tL) H-l/2 

(y + x) 2 



dydx 
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+ H 



x 



i/2-m 



o 



D \ -H+l/2 
X 



t L 



-l ( D V x k ( D \i 

+ x h\ \ — , x ) ax 

\t L J \ \t L J 



D/2t L 



\t L 

D 

x x + | x 

tL 



i/2-m 



where 



h x {D/t L ,x) 



l y H-X,2 {y + D/tL) H-l,2 



y + x 



dy<4/(2H-l). 



Note also that 

.l v ff-l/2 fe + Z , /tl) ff-l/2 



.-/() 



{y + x) 



dy<4x H - 3/2 (D/t L ) H - 



1/2 



so that 

\Z£(t L )\ <4b s t L H C H 
x \{D/t L ) 



H-l/2 



D/2t L 



X 



H-l-rn 



{D/t L -x)- H+1 / 2 dx 



+ (D/t L ) H ~^ 2 f D/tL x- l {D/t L - x) 1/2 -^ dx 

JD/2t L 



(19) 



+ 



+ 



D/2t L 



D/t L 



•}> 2 -^ (D/tL ~ xy H+1 / 2 [x~ l + (D/t L - x)" 1 } dx 



- H+1 > 2 {D/t L - x) l l 2 - r * [x~ l + (D/tL - x)- 1 ] dx\. 

lD/2t L J 

< Ab 5 t L H C H ■ lh(D/t L ) l ~ H ^ 

< Cl tf+H-a-H-r)(l-H- Vl ) 

for an appropriate constant c\. 
Similarly, on A(b), 

\z£(t L )-z£(t L -b 2 -V H s)\ 

= [ D [g'(D, t L , -v) - g'(D, t L - b 2 ^ H s, -v)}Zg(-v) dv 
Jo 

fD/t L 

< I \[g'{D/tL,h-x)-g'(D/tLA-s/t L ,-x)}Zg{-xt L )\dx 
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<b°t L 



5+ H 



/ \g'(D/t L , 1, -x) - g'(D/t L , 1 - s/t L , -x)\x H ~^ dx 
Jo 

f D/t L 

+ \g'(D/t L ,l,-x) 

JD/2t L 

- g'(D/t L , 1 - s/t L , -x)\(D/t L - x) H ~^ dx 



where ii — b 2+1 ^ H ti. Rewriting, 

\zi(t L )-z£(t L -b 2 -V H s)\ 

<b 5 t L H C H (H- 1/2) 



2H- 1 



D/2t L 



X 



1/2-7J! 



D N -H+l/2 
X 



+ 



2H - 1 J D/2 t 



D/ti 



1V(1- S /^) y H - 1 / 2 (y + D/t L ) H - 1 / 2 
lA(l-s/t L ) (y + x) 2 

lV(l-./f £ ) + D/tL) H-l/2 



dydx 



+ 



D/2t L 



.1/2-771 



lA(l-s/ti) 
X 



(y + x) 



dydx 



D 



t X 



h 2 (D/t L ,s/t L ,x)dx 



+ 



D/t L 
D/2t L 



-H+l/2 



l/2-r?i 



D 



x 1 + ( X 

tL 



where 

h 2 (D/t L ,s/t L ,x) 
Since also 



h 2 (D/t L ,s/t L ,x)dx \, 
1V(1-«/*l) y H ~ l / 2 (y + D/t L ) H - 1 / 2 



iA(i-s/ti) y + x 

lV(l--/fi) yH-l/2( y + D / tl )H-l/2 

iiA(i- s /t L ) (y + x) 2 

it follows that 

\m*L)-mL-^ H s)\ 



dy<2\s\/t L . 



dy<2\s\/t L , 
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< 2st7 1 b S t L H C H 



1 L 

| J° /2tL x 1 ' 2 ^ ( D /t L - x)- H+l ' 2 dx 
+ [ D/tL x~ H+l / 2 (D/t L -xf 2 ^dx 

JD/2t L 
rD/2t L 

+ / x l l 2 -^{D/t L -xY H+l / 2 [x^ + {D/t L -xY l }dx 
Jo 

+ [ D/tL x- H+1/2 {D/t L - x) 1/2 ~ Vl [aT 1 + (D/t L - x)- 1 } dx\ 

JD/2t L ) 

<2stl 1 b s t L H C H (D/t L ) 1 - H -^ 

< crfS+H-il-H+rXl-H-rn^l/H-l) 

for an appropriate constant c<i and all large enough b. 

To analyze the factor (021 /an — 1) in the first term on the RHS in (9), 
note that 

ogi = EZ(t L )Z(t L + b 2 ^ H s) - EZJ(t L )ZJ(t L + b 2 - 1 / H s) 
an EZ(t L ) 2 -EZJ(t L ) 2 

The covariance EZ^ (a)Z^ (a + x) satisfies, 

EZ J (a)Z J (a + x) = a 2H f{D/a, 0, x) 

for \x\ < a, where 

f(d, w,x) = H / g(d, 1 + w, -u)[(l + w + x + u) 2H_1 - u 2 ^" 1 ] 
J 

(see Mannersalo [25]). It follows that 

021 ^ _v(8/i L )-l + f(D/t L ,0 i 0)-f(D/t In 0,s/h) 



(20) 



an / 1 -/(£>/*£, 0,0) 

where v(x) = ~(1 + (1 + x) 2H - \x\ 2H ), \x\ < 1. Note that (v(x) - l)/x^ H 
as x — > and it can be easily shown that / (D/tL, 0,0) — /(D/t^, 0, s/^l) = 
o(s/ti,) as 6^ 00, whence there exists a constant C3 such that 

(21) ^-1 <c 3 | S |6-^-D. 

an 

We can then use the above bounds (19)-(21) on the RHS of (9) to conclude 
that there exists C4 such that, on A(b), 

\E(Y b (s)\H) -E(Y h (s)\g)\ < CAb -2H+l b 5+H-(l-H+r)(l-H- m )-(l/H-l) 
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where < e 3 < (1 - H) 2 (l + H)/H - 5 - (1 - H)(r + m ) + rr/i > 0. □ 

Lemma 3. Fix s,t>0. If r < (1 — H) 2 /H then there exists e 2 > such 
that 



\cov(Y b (t),Y b (s)\H) - (6'/2)(\t\ 2 » +\s 
for all large enough b. 



\2H 



\t-s 



2H\ 



< b 



-£2 



Proof. Note from (7) that 
cov (Y b {s),Y b {t)\H) 
= b^- 2H ^ 2 a 32 a 22 

= b 2 ^ 2H ^9 2 EZ(t L + b 2 ~ x l H s)Z{t L + b 2 ~ x l H t) 

- b 2 ^ 2H ^6 2 EZ J (t L + b 2 - l ' H s)Z J {t L + b 2 ~ l l H t) 

_ b 2(i-2H) 6 2 . (EZ(t L )Z(t L + b 2 ~^s) - EZJ(t L )ZJ(t L + b 2 ~ l l H s)) 

varZ(t^) — varZ^(i£,) 

x (EZ(t L )Z(t L + b 2 ~ l ' H t) - EZ J (t L )Z J {t L + b 2 ~ l l H t)) 

= 1^L_((1 + S /t L ) 2 " + (1 + t/t L ) 2H - (|t - S |/t L ) 2 ") 



-^ 2 /(Z?/t L , S /t L ,(t- S )/^) 



? 2 t 2 " 



(«(*/*£) - f(D/t L As/i L ))(v(t/t L ) ~ f(D/t L ,0,t/t L )) 



1- f(D/t L ,0,0) 

where v (•) and /(•, •, •) are as above. Straightforward algebraic manipulations 
then give 

[cov(f fe ( S ), Y b (t)\H) - cov(Y(s), Y(t))]6~ 2 tZ 2H (l - f(D/t L ,0,0)) 
1-/(5/^,0,0) 



1 + 

t L 



211 



2H 



+ (l +t /t L ) 2H -\t/t L \ 2H -2f 
( D 



(D s . 



f 



,0,^- 

V*L t L 



v(t/t L )-f( —,0,t/t L 



D 



t L 



(1-/(5/^,0,0)) 



f(D/t L ,0,t/h) + f[D/t L ,0, 



t L 
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f(D/t L , 0, 0) - / ( D/t L , j-,(t- s)/t L 



f[D/t L ,0,j-)-f(D/t L ,0,0) 



[f(D/t Ll 0,t/h)-f(D/t L ,0,l 



(v(t/i L ) - 1) 



t L 



f[D/t L ,0,^)-f(D/t L ,0,0) 



l)(v(t/t L )-l), 



and the last four terms can be shown to be 0(b 2< ^/ H 1 )) = 0(i L 2 ) as b- 
oo, so that 

[cov(Y\s),Y b (t)\H)-cov(Y(s),Y(t))]e- 2 tl 2H 

= f(D/t L ,0,t/t L )+f(D/t L ,0,s/t L ) ~ f(D/t L ,0,0) 
- f(D/t L ,s/t L , (t - 8)/i L ) + 0{tl 2 ) 

"D/tL 



Hi {g(D/t L ,l,-u)[(l + u + t/t L ) 
lo 



2H-1 



+ (l + u + s/t L 



2H-l +u 2H-l 



g(D/t L ,l,-u)[(l + u) 



- g(D/t L , 1 + s/t L , -u)[(l + u + t/t L ) 

+ o(i- L 2 ) 



2H-1 „.2H-1 



u ZH ~ l ]}du 



H 



D/t L 



{g(D/t L , 1, -u)[(l + u + S /tL? H - 1 - (1 + u) 2H ~ l ] 
-(l + u + t/h) 211 - 1 



x [g(D/t L , 1 + s/t L , -u) - g(D/t L , 1, -u)]}du 
+ 0(t- L 1 Dtl 1 ), 

where the last step follows because 



D/t L 



2H-1 



g(D/t L , 1, -u) - g(D/t L , 1 + s/t L , -u)\ du 



rD/t L 

C H / u H ~^ 2 (D/tL — u)~ 
Jo 



H+l/2 
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-lV(l+«/t £ ) yH-l/2( y + D / tL }H-l/2 



lA(l+s/ti) 



y + u 



dydu 



[ D ' tL u H-l/2 {D/tL _ u) -H + l/2 du 
JO 



<2\s/t L 

= 2\s/t L \{D/t L 
= 0{tl 1 Dtl 1 ) 



.H- 



as b — ► oo. Put 



h(d, x, u) 



l+x 



y H-l/2 {y + d) H-l/2 



dy, 



y + u 

and note that there exist constants cq and cj such that, for all small enough 
d and all < u, \x\ < d, 

\(l + u + x) 2H - 1 - (1 + u) 211 - 1 - {2H - 1)(1 + u) 2H ~ 2 x\ < c 6 x 2 , 

\h(d,x,u) - h(d,0,u) - x(l + d) H ~ 1/2 (l + u^l < c 7 x 2 . 

It follows that, for all large enough b, 

\cov{Y b {s),Y b (t)\H)-cov(Y{s),Y{t))\6- 2 t7 2H 



H 



D/t L 



U 



-H+l/l( 



D 



\t L 



-H+l/2 



[(1 +u + s/h) 211 - 1 - (1 + u) 2 "- 1 ^ 0, u) 

_(l +u + t /t £ )2H-l 

x [h(D/t L ,s/t L , -u) - h(D/t L ,0, -u)} 



du 



rD/t L 

H ' 



ii 



H+l/2 {D/tL _ u) -H + l/2 

x\(2H- l)(l + u) 2H - 2 h(D/t L ,0,u)(s/t L ) 

_ (i + u + t/t L ) 2H -\l + D/t L ) H - l / 2 {\ + u)-\s/t L )\ du 
+ 0((s/t L ) 2 (D/t L ) 2 ^) + O&Dti 1 ). 

But there exist constants cs, eg such that, for all small enough d and all < u, 

\y\ <d, 



(l+d) H ~ 1/2 -l<c 8 d and Kl + n + y) 2 ^" 1 - (!+■ 



\2H-1\ 



< cg\y\, 
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so that, for all large enough b, 
\cov(Y b (s),Y b (t)\H)-cov(Y(s),Y(t))\e- 2 tl 2H 



< H{s/i L ) 



D/t \-H + l/2fD 



-H+l/2 



(2tf-l)^£,0,u) -1 



(1 + u) 



2H-2 



du 



+ 0((s/t L )(D/t L f- 2H ) + 0((s/t L ) 2 (D/t L ) 2 ^) + 0(r L l Dt- L l ). 

But, for all small enough d > 0, \{2H - l)h(d, 0, u) - 1| < 5(1 - H)~ l d 2H ~ l , 
so that for an appropriate constant cio and all large enough b, 

\cov(Y\s),Y b (t)\H)-cov(Y(s),YmO-% 2H 

■ D \-H+l/2 



< ci (s/t L ) 



n \2H-l r D/t L 

TJ ' u 



(l + u) 



2H-2 



du 



+ 0(t^Dt L > 

= c 1Q {s/t L ){D/t L f H - l O{{D/t L f^- H )) + Oit^Dtl 1 ) 
= 0(r L 1 (D/t L )), 
and we conclude that 

\ C ov(Y b (s),Y b (t)\n)-cov(Y(s),Y(t))\=0(P L H ~HD/t L )) 

= 0{b ^l~Hf/H + r ) 

= o{b- £2 ) 

for any < e 2 < -{I - H) 2 /H + r. □ 

Lemma 4. Assume (5) and (4) hold. Given M > 0, for all large enough 
b the function E{Y b {u)\TL) is differentiate on [—M,M] and there exists 
a constant Cm such that, on A(b), 

dE{Y b (u)\H) 



(22) 
ue [— M, M]. 



du 



< Cm, 



Proof. Recall that, with / and v as in the proof of Lemma 3, 
•rhi \\<u\ _ ul-2H 1^/^,0, u/t L ) - v(u/t L ) 



E{Y b {u)\H) = b l 



1 0,0) 



(23) 



+ b 1 - 2H z£(t L + b 2 - 1 / H u) 

+ [E{Y b (u)\G) - E(Y b (u)\IC)] + E(Y b (u)\IC). 
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The first term on the RHS in (23) is differentiable (as a function of u), and 
d /^ 2H+1 f(D/t L ,0,u/t L )-v(u/t L ) 



du 



li 



l-/(D/i L ,0,0) Q 
(df/dx 3 )(D/t L ,0,u/t L ) - v'(u/t L 



l1-2H 



l-f(D/t L ,0,0) 

Here, the first factor on the RHS is bounded over [—M,M] for large enough 
b, while the second factor does not depend on u and is o(l) as b — > oo, as can 
be seen from (19) in the proof of Lemma 2. Hence, there exists a constant 
Cmi such that, on A(b), 



ou 



v(u/t L )- f(D/t L ,0,u/t L ) ~ Ju 



< Cmi, 



l-f(D/t L ,0,0) 
ue[-M,M]. 

For the second term on the RHS in (23), the integral representation used 
in the proof of Lemma 2 can be used to obtain 

d 



A-2H 



= b 



Z$(t L + b^/ H i 



hCh ( D/tL „i/2~h(D 

tL 



ti Jo 



X 



1/2-H 



(1 + u/t L ) H - 1/2 (l + D/t L + u/t L ) H - 1/2 

(l + x + u/i L ) 
H-l/2 H-l/2 1 



Zg{xti) dx. 



X D/tL — X 1 + x + u/ti. 

Using that, on A(b), \Zg(-w)\ < max(w H b 5 ,w H ~ r ' 1 ) for < w < D/2 and 
\Z g (-w)\ <max((D -w) H b 5 ,(D -w) 11 - 111 ) for D/2<w<D, one can con- 
clude that the RHS on the above equation is bounded, so that there exists 
Cm2 such that, on A(b), 

b^ H \^{t L -b^l H u)\<C M 2, 

ue [—M,M]. 

For the third term on the RHS in (23), recall 

E(Y\u)\G)-E(Y b (u)\lC) 

1-2H Z K {-t R ) 



9b 



4(1 + D/t L y H - (1 + (1 + D/t L f H - {D/t L f H f 

2H 







u 


2H 




\ 2H ( D 


u\ 


X 


2 




+ 2| 




) -i i+ r L 








T L 
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+ 1 



tj 



2H 



D 

l + r L 



211 



If) 



whence 



— (E(Y b (u)\g)-E(Y b (u)\IC)) 



2Hf L i 6b 



laa-2H 



4(1 + D/t L y H - (i + (i + D/t L y H - {D/t L y H y 



2sign(u) 



/ D u_^ 2 

tL t L , 



D ^2H 

V + T L , 



t L 



x I sign(u) 



2H-1 



1 + 



u 



2H-1 



< {l + M 2H ~ l )0{b- £2 ), 



where the term 0(b £2 ) does not depend on u, and £2 is as in Step 2. Hence, 
there exists a constant Cm3 such that, on A(b), 

d 



du 



(E(Y\u)\g)-E(Y b (u)\!C)) 



< Cm3j 



ue[-M,M}. 

For the last term on the RHS in (23), one can differentiate the expression 
for E(Y b (u)\)C) given in (3) (with u in place of s) to obtain 

^-E{Y b {u)\K) = -2^K 1 |n| 2H - 1 sign( U )(t7t L ) 2H 
ou 



2-l h l-H+r 



, + Utr b 

(l-H)t L y y ' 

+ Hb l ~ 2H tl l [X{-t L ) + b + fib H+r ] 

x (lu/til^-isigii^) - (1 + u/i L ) 2H - 1 ), 

and since \X(— ix) + b\ < W < b H+r on A(b), we conclude there exists a 
constant Cma such that, on A(b), 



^-E{Y b {u)\K) 



< Cm4, 



u G [—M,M]. Taking Cm = Cmi + Cm2 + Cm3 + CWi the result follows. 

□ 



28 H. AWAD AND P. GLYNN 

Lemma 5. i/7, Tji, r and 5 satisfy (15)-(17), then 
P{A(b)A{T b < 00}} = o(P{T b < 00}). 

PROOF. We show that all the terms on the RHS in (14) are o(P(T b < 
00)) as b — ► 00. 

For the first term, note P(A 5 ) < P(T b < t L ) = o(P(T b < 00)) by Theo- 
rem 3. 

For the second term, choose any e such that < e < (2H — 1)(1 — H — r)/8, 
and put vi = 1 - H + if 2 + rH + e. Note that 

p(A)n{XH L )>-&+&" 1 })<p( sup z{-t L )-z{-t L -s)>b^/e). 

and it then follows from Borell's inequality (Lemma 7 below) that there 
exist constants c\ and C2 such that 

P{A n {X(-t L ) >-b + b U1 }) < a exp(-c 2 6 2iyi /b 2H{ - H+r ^) = o{P{T b < 00)), 

where the last follows since 2v\ — 2H(r + H) > 2(1 — if) and 

P(T b < 00) ~ g(6) exp(-(«2/(2(l - H)))b 2 ^ H ^), 

where q(-) is a ratio of polynomials, by a result of Hiisler and Piterbarg [15], 
Corollary 2. 

Let v\ and e be as above, and for k > 1 put 

= (if + r - (2H - 1)(1 - if - r))/2 + (y k V (if + r))/2 + e. 

Observe that fi > r + if, and it follows from Lemma 6 that 

P(A) n (X(-ti) > -b + b^}) = (P(T b < 00)) 

as b— > 00. Noting that u k is nonincreasing in k, an inductive application of 
Lemma 6 shows that whenever v k > H + r, then 

P(A n {X{-t L ) >-b + &•*+!}) = o(P(T b < 00)). 

Also, note that if k > 2 and > if + r, then 

i/fc - i/fe+i = (vk-i ~ v k )l2. 

It then follows by an inductive argument that, if Vf, > if + r, then 

z/ fe+1 =H + r-(2H - 1)(1 - if - r) + 2~ fc if (1 - if - r) + 2(1 - 2" fc+1 )e. 

Let fe* = inf{A;:i/jt < ff + r} and note 1 < k* < 00. Also, note that ffc* > 
if + r-(2if-l)(l-if-r)/2. Then, putting v* = v k * + i = H + r - (H - 
1/2) (1 — if — r) + e, it follows by Lemma 6 that 

P(A) n {X(-i L ) > -& + = o(P(T & < 00)). 
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If 7 satisfies (15) and e is small enough, then v* < 7, so that the second 
term on the RHS in (14) satisfies 

P(A) n A*) = P(A n {X(-t L ) >-b + W}) 

< P(A n {x(-t L ) >-b + b v '}) = o(P(T b < 00)). 
For the third term on the RHS in (14), note that 

P(Ain^) 

= EI Al P(\Z K (-t R )\>b^\IC) 

< P(X(-t L ) <-b + V) ■ P(2 H b H ^ H+r ^\N(0, 1)| > ^ 3 ) 

But it follows from the Hiisler-Piterbarg asymptotic (13) and equation (26) 
in the Proof of Lemma 6 that 

P(X(-t L ) <-b + 6 7 ) = P{T b < oo)q(b) exp{b 2(1 - H) 0(b- 1+H+r )), 

while 

P(2 H b H ^ H+r ~>\N(0,l)\ > b^) = 0(b^ +H ^ H+r hM-b 2{m ~ H2 ~ rH) /2 2H - 1 )) 

in view of (16), it then follows that P(A 1 n A§) = o(P(T b < 00)). 

To deal with the fourth term on the RHS in (14), put TV = \W e ], h = 
b -2(i-H)/ Vl , an d A = (t* - t L - h)/N, where % > 0. Note that 

N N 

p(A x r\A$)<Y, p (^i n E n ) + ]T P(Ai n B n ) 

, n=0 n=0 

24 

00 00 

+ P(^i n c n ) + p (^i n °n), 

71=1 71=2 

where 

E n 4 {\Zg(-t L -h- nA)\ > (l/2)b S (h + nA) H }, 



= { sup \Zg(-t L -h-nA)-Zg(-t L -h-nA-s)\>(l/2)b 5 h H \, 
C n 4 {\Zg(-t L - h/n 2 )\ > (l/2)(h/n 2 ) H ~ V1 }, 

D n ±{ sup \Zg(-t L -s)-Zg(-t L -h/n 2 )\>(l/2)(h/n 2 ) H -^}. 

l hn~ 2 <s<h{n-l)- 2 ' 

Note 

P{A X n £ n ) 

= E/ Al E(/ £ J£) 
(25) < EI Al P((h + nA) H |iV(0, 1)| > (l/2)(/i + nA) H b s ) 

< P(X(-t L ) <-b + b^)0(b- 5 exp(-6 25 /8)) 
= P(T 6 < oo)g(6) exp(6 2 ( 1 - i/) 0(6- 1+H+r ))0(6- 5 exp(-6 2 V8)), 
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where q(b) is a ratio of polynomials and the last step follows from the Hiisler- 
Piterbarg asymptotic (13) and equation (26) in the Proof of Lemma 6. Since 
S > (1 — H + r)/2, it follows that there exists a constant c 3 > independent 
of n such that 

P(Ai n E n ) < P(T b < oo)g(6)6-' 5 0(exp(0(6 1 ~ //+r ) - b 28 /8)), 
<P(T 6 <oo)exp(-c 3 6 25 ). 

Hence, 

N 

P(4i n E n ) < P(T b < oo)2W« exp(-c 3 & 25 ) = o(P(T b < 00)) 

n=0 

as b — ► 00. 

Next, note that on A\, 

\Z G {-t L -h-nA + s)- Z G (-t L -h-nA)\< c 4 A 

for some constant C4 > 0, all < s < A and all < n < N. Thus, provided 
7/6 is large enough so that = o(b s h H ) as b — > 00, it follows that for large 
enough b 

PiAtHBn) 

< P( sup \Zg(-t L -h-nA-s)- Z g (-t L -h-nA)\> (l/A)b s h H ) 

<c 5 eM-c G b 25 h 2H A~ 2H ), 

for appropriate positive constants C5 and cq, by Borell's inequality. Hence, 
by choosing t/q large enough we can ensure 

P(i 1 nB n )<exp(-6 4 < 1 ^) ) 

and it then follows that 

N 

P(^i n B n ) < 2W 6 exp(-6 4 ( 1 ~ H )) = o(P(T b < 00)) 

n=0 

as b — ► 00. 

A similar argument as the one used to analyze ~P(A\ E n ) shows that 
there exist positive constants ci, 02,03 such that 

P ( Ax n C n ) < ci exp ( - c 2 /i ~ 2r?1 ) exp ( - c 3 ?i 4??1 ) . 

Also, a similar argument as the one used to analyze P(^4i H B n ) shows that 
there exist positive constants 04,05, cq such that 

P(Ai n D n ) < 04 exp(-c 5 /i" 2m ) exp(-c 6 n 4r ' 1 ). 
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It then follows that the last two summations in (24) are of order 0(exp(— (c2 A 
c 5 )h- 2 ^))=o(P(T b < oo)), so that P (A 1 n ) = o(P (T b < oo)) as b^oo. 
The proof that P(AiHA^) = o(P(T b < oo)) uses the same type of argument, 
and is omitted for brevity. 

Hence, all terms on the RHS in (14) are o(P(T b < oo)), what had to be 
shown. □ 

Lemma 6. Let r be the small number chosen in the definition of Aq. If 
1 > v > H + r - (2iT - 1)(1 - H - r)/2 and P(A n {X(-t L ) >-b + b u }) = 
o(P(T b < oo)), then for any v > (H + r- (2H- 1)(1 - H - r))/2+ (i/V (H + 
r))/2 it also holds that 

P(A (b) n {X(-t L ) >-b + b°}) = o(P(T b < oo)) 

as 6 — > oo . 

PROOF. Observe that 

p(A (b)n{x(-t L )>-b + b>}) 

= P(A (b) n {X(-t L ) G [-b + b°, -b + b u ]}) + o(P(T b < oo)) 
= E[I(X(-t L ) G l-b + b*,b- b v \) ■ P(A (b)\K)} + o(P(T b < oo)). 

But on {X{-t L ) G [-b + b°,-b + b u ]}, E(X(-t L ) - X(-t L - s)\lC) = o(l) as 
b— > oo. In particular, for all large enough b, 

I(X(-t L )e[-b + b*,-b + b»])-l Ao{b) 

<l( sup -Z K (-t L -s)>b»/(28)). 

By Borell's inequality, we can find positive constants c\ and C2 such that 
P( sup -Z K {-t L -s)>U'/{2e)\K) < Cl exp(-c 2 6 2 76 2 ^ +r )) 

on {X(t L ) G [-b + tf ,-b + b u ]}. It then follows that 
P(A (b)n{X(-t L )>-b + b>}) 

< P(X(-t L ) <-b + b u ) ■ a eM-C2b 2 *~ 2H{H+r) ) + o(P(T b < oo)). 

But 

P(X(-t L )<-b + b u ) 

(26) 

~ 9(1 - J f/)6- 1 exp|-^ I ^6 2 ( 1 ^)(l - *»(&))}, 
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where 

a 1 _ {l-(^/{l-H))b H +^-b^/(l-H)f _ 1+uV(H+r) 
as b — > oo. Hence, with q(-) a ratio of polynomials, 

p(A (b)n{x(-t L )>-b + b i> }) 

< P(T b < oo){g(b) ex p( 2 ^ ( ff ) &2(1 " H) ) exp(-c 2 6 2 ' > - 2 ^+ r )) + o(l)} 

< P(T 6 < oo){g(&) exp(0(6 1 " 2H+w ^ +r )) - C2b 29 ~ 2H ^ H+r ^) + o(l)} 
= o(P(T 6 <oo)) 

as b -> oo, since 1 - 2iJ + z/ V (il + r) < 2P - 2H(H + r). □ 

The following result is known as Borell's inequality in the literature; see 
Adler [1], Theorem 2.1. We quote it here since we use it at many points in 
the proofs of Theorems 1 and 3. 

Lemma 7 (Borell's inequality). Let {V(t):t£T} be a centered Gaus- 
sian process with sample paths bounded a.s. Let V M = sup tg y V(t). Then 
EV M < oo and for all A > 

P(\V M - EV M \ > A) < 2e -( 1 /2)A 2 /4, 

where a\ = sup tg y var V(t). As an immediate consequence, for A > ~EV M , 

P(V M > A) < 2 e -( 1 /2)(A-EV«)V4. 

In particular, if Z is standard fBM and C = Esupg <s<1 Z(s), then for 
\>t H C 

Pf sup Z( S )>A)<2e^2)(A^cW". 

Table 1 summarizes selected notation which is used repeatedly at different 
points of the argument. 
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tt = t*-b H+r ; t R = t*+b H+r ; t L = b- 2+1 ? H t L ; D = 2b H+r 

Y\t) =b- 2H+1 X(b 2 - 1 / H t) 

Z(t) =Z(-t L +t)-Z(-U) 

K =a{Z(-t L ))=a(Z{t L )) 

g = a {Z(-t L ),Z(-t R )) = o-{Z{t L )^Z{-D)) 

H = er(Z(a) : - t R < s < -t L ) = o(Z(t L ), (Z{s) : - D< s < 0)) 

J =a(Z(s): -D<s<0) 

Ao(b) = {3« G l-t R , -t L ] : X(s) = -6} 

Ai(b) = {X(-t L )€[-b,-b + b' 1 ]} 

A 2 {b) = {\Z K {-t R )\<W 3 } 

A 3 (b) = {\Z g (-t L -s)\< max(s H 6 i ,s fl -" 1 ),0< s < b H+r } 
A A [b) = {\Z g \-t R + s)\< ma,x(s H b s , s H -^),0< s < b H+r } 

A{b) = A (b) n Ai(b) n A 2 (b) n A 3 (b) n A 4 (b) 
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